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A role of umklapp scattering has been examined for two-coupled chains with both forward and backward 
scatterings by applying renormalization group method to bosonized Hamiltonian. It has been found that a 
state with relevant interchain hopping changes into a state with irrelevant (confined) one when the magnitude 
of umklapp scattering becomes larger than that of interchain hopping. Critical value of umklapp scattering 
for such a confinement-deconfinement transition is calculated as the function of interchain hopping and 
intrachain interactions. A crossover from one-dimensional regime into that of coupled chains is also shown 
with decreasing temperature. 

PACS numbers: 71.10.Hf, 71.10.Pm, 71.30.+h, 75.30. Fv 



I. INTRODUCTION 

Quasi-one-dimensional organic conductors, (TMTSF^X 
and (TMTTF)2X salts, exhibit instabilities toward spin 
Peierls state, spin density wave (SDW) state and supercon- 
ducting (SC) state, where the phase diagram has been dis-. 
played on the plane of effective pressures and temperatureoEl 
The interplay of low-dimensionality and repulsive inter- 
action is important for the SDW state which indicates one- 
dimensional fluctuations .Era There are also some evidences for 
dimensional crossover^ 

Crystal structure shows quarter-filling for conduction elec- 
trons rput the existence of dimerization leads to half-filled 
band.B A crossover from half-filled band to quarter-filled one 
has been found by decreasing dimerization under effective 
pressure, i.e., the variation of anions, X. Electronic properties, 
which suggest a role of the dimerization, have been reported 
recently at temperatures just above the SDW state.Q'Q Op- 
tical experiments on a series of above materials, which have 
different values of interchain electron transfer energy, show a 
correlation gap due to umklapp scattering and a crossover 
from metallic state to insulating state with increasing the 
anisotropy. An insulator to metal transition followed by the 
deconfinement of interchain hopping has been observed when 
the interchain transfer energy exceeds a critical value with a 
magnitude of the order of the gap. 

Theoretical studies of these conductors have been explored 
by use of quasi-one-dimensional model consisting of an ar- 
ray of chains coupled by interchain hopping. For repulsive 
intrachain interaction and incommensurate band, the trans- 
verse hopping is always relevant for the weak interactions but 
there is a reduction of transverse hopping by one-dimensional 
fluctuations Two-coupled chains is a basic model for quasi- 
one-dimensional system since both one-dimensional fluctu- 
ation and transverse hopping can be studied on the same 
footing. In a Tomonaga-Luttinger model with only forward 
scattering, the dominant state remains the same as that of 
one- dimensional system but the degeneracy of in-phase and 
out-of-phase pairings is removed.LJLJ When backward scat- 
tering is added, the phase diagram becomes quite different 
from that of a single chain. In a Hubbard model with repul- 
sive interaction and without umklapp scattering, the ground 



state of two-coupled chains is the d-wave like SC stateE10~ll3 
although that of a single chain is the SDW state.EJ The ef- 
fect of interchain hopping is much strong compared with the 
intrachain interaction since the transverse hopping-Ls relevant 
except for extremely large intrachain interaction.E3 However 
intrachain interaction becomes important as well as the inter- 
chain hopping for the case of the spin anisotropic backward 
scattering where a spin gap induced in a single chain leads to 
a competition between SDW state and SC state£j 

Recently confinement which denotes incoherence of sin- 
gle particle, hopping between Luttinger liquids has been 
maintained J2JEj where there is no coherence of hopping and 
then no split Fermi surface below a critical value of single par- 
ticle hopping. The confinement has been argued for the metal- 
lic state of organic conductor (TMTSF^X under a magJietic 
field, which is close to coherence-incoherence transitionEd The 
role of umklapp scattering, which leads to the relevance and 
the irrelevance of the pGprrelation gap, has been examined 
for organic conductors H3 In terms of a Mott gap, the irrele- 
vance of single particle bowsing has been discussed in quasi- 
one- dimensional systemBi'LJ A confinement has beea demon- 
strated in two-coupled chains with half-filled bandtj in order 
to understand a crossover from metallic state to insulating 
state, which has been found at temperatures just above SDW 
state of organic conductors.^ The weakly coupled half-filled 
chains with infinite numbers have been also-studied by a per- 
turbative renormalization group approachE3'E3 

In the present paper, such a deconfinement-confinement 
transition due to umklapp scattering is studied in details for 
two-coupled chains with half-filled band by developing the 
previous work.E3 In section II, formulation is given in terms of 
bosonized phase Hamiltonian. Renormalization group equa- 
tions are derived for coupling constants and response func- 
tions. In section III, the critical value for confinement is cal- 
culated. A crossover at finite temperatures is also shown. In 
section IV, we discuss the validity of our present calculation 
and examine an effect of forward scattering within the same 
branch. 

II. FORMULATION 



We consider two-coupled chains given by 
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where t is the interchain hoping energy. The quantity a\ a i 
denotes a creation operator for the electron with momentum 
k, spin cr(=t, J. or + , — ) and chain index i(= 1, 2). The sym- 
bol p _ +(— ) represents the right-going (left-going) state. In 
Eq. (2.1), tk.p(= v?(pk — kp)) is the linearized kinetic energy 
with Fermi velocity vf and Fermi momentum kp. Quantities 
g 2 , gi and g$ are coupling constants of intrachain interactions 
for forward scattering, backward scattering and umklapp scat- 
tering, respectively. 

The diagonalization of the first and the second terms in 
Eq. ( |2.l| ) is performed by making use of an unitary transfor- 
mation, Ck, p ,<T,p = (— fJ.a k ,p,<T,i + afc, P , CT ,2)/\/2 with fj, = ±. Af- 
ter the bosonization of electrons around the new Fermi point, 
= kp — fit/vp, we define the phase variables, 9 P + and 9 a + 
(9c+ and 6s+), which express fluctuations of the total (trans- 
verse) charge density and spin density.E-3 They are given by 
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There is a commutation relation that [0 V + (x), 9 v i- (a/)] = 
i^&v.v 1 sgn(x — x) where the suffix — denotes the canonically 
conjugate variable. In terms of these phase variables, the field 
operator is expressed as 
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where a is of the orde r of the lattice constant. The phase 
factor, 7rH Pj(TiM , in Eq. ( |2.6| ), which is introduced for the anti- 
commutation relation, is taken as 

(-1V+ 1 . 

E 2n+j =N!+---+N 2n + - — ^ (N 2 n+i + N 2n+2 ), 

(2.8) 

where j'=l,2 and n — 0,1,2,3. The quantity, Ni, denotes 
number operator, and the suffix i is related to (p, a, /i) 
as, (+,+,+) = 1, (+,-,+) = 2, (+,+,-) = 3, (+,-,-) = 
4, (-, +, +) = 5, (-, -, +) = 6, (-, +, -) = 7 and (-, -) = 
8, respectively. Such a choice of H Pj(TiM conserves a sign of 
interactions, which are repres ente d by phase operators. In 
terms of these operators, Eq. (2.1) is rewritten as (Appendix 
A) 
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where = vpy/l — {g v /2irVY} 2 , Kv = [{1 — g» /2tvvf} / '{1 + 
g^Ti-vp} 1/2 , p p = 2g 2 - gi, g a = -gi and g c = gs = 0. 

We maka i*se of renormalization group method for response 
functionsJ^jEj which are assumed to be invariant for scal- 
ing, a — *■ a' — ae . We express the nonlinear terms in 
Eq. (2.£) as g vp y p i/{2ir 2 a 2 ) J dx cos ^/29 vp cos \f28 v i p i where 
\fW U p = V%6vp — Atx/vF for v = C and p = + and 
v20up = \/26 vp otherwise. Then the coupling constants 
are given by gc+,s+ = -gc-,s- = (g P + fla)/2, ffc+,s- = 
-5c-,s+ = (ff P - g a )/2, g a +,c+ = -fftr+,c- = -ff<x+,s+ = 
-5<r+,s- = ffi and g p+ ,c+ = ff P +,c- = -g P +,s+ = 5p+,s- = 
g3. Response functions defined by Ra(xi — %2,ti — T2) = 
(TtOa (an, ri)0^ (* a , r 2 )) are evaluated for SDW and SC 
states where Tj is the imaginary time and Oa denotes the 
order parameter. Then renormalization group equations are 
expressed asf— I (Appendix B) 
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where t(l) = t(t)/ep, ef = hfQ -1 , Si/ = v„/vf, v — p, 
a and p, p' = ±. In these equations, the /-dependence is 
not written explicitly and J„ (n = 0, 1) is the n-th order 
Bessel function." Initial conditions are given by K v (0) = K v , 
G uPlV i p i(Q) = g„ PtV , p r /2ttv f and t(Q) = t/e F . 

The second order renormalization group equations with re- 
spect to all the coupling constants are derivedrbjy expanding 
as K^ 1 ^) — 1 =p G V {1) + • • •. In case of gs-f= 0,t3 these equa- 
tions become equal to those of Fabriziojia which satisfy the 
SU(2) symmetry with respect to spin rotation. Althoug h suc h 
a, sym metry is satisfied only approximately for Eqs. ( |2.10| )- 
(I2T2J), the difference is very small within the present choice 
of parameters as is s hown later. The renormalization equa- 
tions of Eqs. ( 2. 10 )-( 2.12 ) determine the fluctuations of the 



total charge, total spin, transve rse ch arge and transverse spin 
density respectively. Equation (2.15) corresponds to forward 
scattering and backward scattering with parallel spins. In 
r.h.s. of these equations, there are bilinear terms with re- 
spect to G vp>u ipi(l), which appear in the presence of umk- 
lapp scattering and/or backward scattering whil e the y are 
absen t for only forward scattering. Equations (2.13) and 
( |2.14[ ) with v = p (y = a) correspond to umklapp scattering 
(backward scattering with opposite spins). It is found that 
there is a symmetry between equations of the total charge 
and those of the total spin, i.e., the renormalization equa- 
tions remain the same for the replacement given by K p ^+ 
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(G CT+ ,c+, G CT+ ,c-, G CT+ ,s+, Gv+,3-). Equation ( |2.16[ ) is the 
scaling equation for the interchain hopping. It is noted that 
these [equations with t = is reduced to those of a single 
chainB 

We examine order parameters for the possible states in case 
of repulsive interaction. In terms of phase variables, order pa- 
rameters are expressed as (Appendix A) 
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where Vwfc) = (1/VZ) E fc e ite a fe>p , CT ,j. In Eqs. (Q- 
( pJ9| ), LSDW||, out (TSDW|| j0ut ) denotes longitudinal (trans- 
verse) SDW with intrachain and out-of-phase pairing. The 
suffix, SSx.in, represents SCd state, i.e., the singlet SC state 
with interchain and in-phase pairing. 

The renormalization group technique is also applied to the 
calcul ation of re sponse functions for the order parameters, 
Eqs. ( 2. 17 )-( 2.1£ ) . Normalized response functions are derived 
as (Appendix B) 
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where r = [x 2 + (uft) 2 ] 1 ' 2 and the quantities K v (l) (u — p, 
a, C and S ) and G uPtV > p r(l) (p, p' — ±) are calculated from 
Eqs. ( 2.1C )-( In these equations, the renormalization 
of the velocitjET has been discarded in a way similar to the 
spinless caseEj 



III. CONFINEMENT VS. DECONFINEMENT 

We examine confinement-deconfinement transition by cal- 
culating the renormalization group equations for interactions 
of both Hubbard model and general model with gi ^ g2- The 
scaling quantity I (= lnr/a) is related to energy ui and/or 
temperature T by the relation that / = \n(eF /lo) = \n(ep/T). 
Numerical calculation is performed by use of normalized 
quantities §j = gj/(2-Kv-p) for j = 1 ~ 3. In Fig. 1(a), quan- 
tities t(l) and 1/Kc(l) as a function of I are shown by solid 
curves and dashed curves, respectively with the fixed §3 = 
0.05, <?3c (=0.119) and 0.3 where t/e F = 0.1 and §1 = §2 = 0.4. 
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Both mutual interactions and umklapp scat tering suppress 
the increase of t(l) as is seen from Eq. ( 2.K ). In case of §3 
=0.05, t(l) (solid curve (1)) increases rapidly. Such a behav- 
ior of t(l) denotes the deconfinement of the transverse hop- 
ping. The corresponding l/Kc(l) shown by dashed curve (I) 
decreases monotonically to zero indicating a formation of the 
transverse charge gap. In the present case, some of G V p v t p t(l) 
diverge at finite / and then solutions stop due to the second 
order renormalization group equations. It is expected that 
the calculation with third ordipi equations gives finite value 
of G uPjI/ i p i(l) for all value of Lt3 A noticeable difference ap- 
pears for large value of the umklapp scattering as is shown for 
i?3 = 0.3 (curves (3) and (III)). With increasing I, t(l) (curve 
(3)) takes a maximum and reduces to zero and 1/Kc(l) (curve 
(III)) remains finite even at the limiting value of I. Such a 
behavior of t(l) indicates the absence of interchain hopping 
which leads to confinement of electrons within a single chain. 
There is no transverse charge gap due to finite Kc(l), where 
the oscillatory behavior comes from the Bessel functions in 
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FIG. 1. (a) The /-dependences of t(l) (solid curves) 
and 1/Kc(l) (dashed curves) for §3= 0.05 ((1) and (I)), 
53 = g-ic(= 0.119) ((2) and (II)) and g 3 = 0.3 ((3) and (III)), 
respectively where t/ep = 0.1 and g\ = §2 = 0.4. In the 
inset, curves (1), (2) and (3) denote t(l)/te l for c/3 = 0.05 
(1), g 3c (2) and 0.3 (3). (b) The ^-dependences of ¥^ and 
1/K^ ynl . The quantity t min denotes a minimum of t(l)/te 
and the quantity 1/A"^f ym is the limiting value of 1/Kc(l)- 
The arrow denotes the critical value, §3 = g 3c , corresponding 
to a boundary between deconfinement and confinement. 



Eqs. ( 2.10 )-( 2. 15 ) obtained with use of the sharp cutoff in 
the formulation of renormalization group technique.Ej The 
fact that G p+ ,c+(0/G P +,c-(0 ^ G (T+ ,c + (0/G CT +,c-(0 ^ 
Go+,s+(0/Go-,s+(Q ^ G c +,s-(0/Gc-,s-(0 =a l/JCo(I) for 
the limiting value is consistent with the irrelevance of the 
interchain hopping. At a critical value given by §3 = §3 C , 
a transition from deconfinement to confinement takes place 
where both t(l) and 1/Kc(l) reduce to zero at the limiting 
value of I (curves (2) and (II)). In the inset, the normalized 
interchain hopping, t(l)/te , is shown for §3 = 0.1 (1), 0.119 



(2) and 0.3(3) where te l denotes the value for non-interacting 
one. The limiting behavior of curve (1), which remains con- 
stant for large I, indicates deconfinement. In Fig. 1(b), the (in- 
dependences of l/K^ 7 " 1 and t min are shown where 1/K^ yl!n 
is the limiting value of Kc- The quantity t mln denotes a min- 
imum of t(l)/te , which is found with increasing I from zero 
(e.g. curve (1) in the inset of Fig. 1(a)), and is essentially 
the same as the limiting value. Deconfinement is obtained 
for finite t mln while confinement is found for finite l/KQ yul . 
Both t min and 1/A"^f ym are reduced to zero at 53 = g3 C , which 
denotes a critical value for deconfinement-confinement transi- 



tion. We note that the Bessel function Ji(4t(£)) in Eq. (2.16) 
is crucial to obtain such a transition. Actually, in r.h.s. of 



Eq. (2.16), the second term becomes negligible for deconfine- 
ment but the second term becomes larger than the first term 
for confinement. 

In Fig. 2, the corresponding /-dependences for K p (l), K a (l) 
and K g (I) are shown by solid curves, dotted curves and dashed 
curves respectively where numerical results are shown for 
\G vp y p >(l)\ < 10. Curves (1), (4) and (7) are for g 3 = 0.05, 
curves (2), (5) and (8) are for §3 = (?3 C and curves (3), (6) 
and (9) are for §3 — 0.3. The quantity K p (l) as a function of 
/ decreases to zero. A charge gap is formed for K p (l) ~ K p /2, 
which [gives a result consistent with that of the Hubbard 
model.LJ The transverse spin fluctuation is also suppressed by 
umklapp scattering because Ks(l) with the fixed / is reduced 
by §3. However the ^-dependence of K a (l) is very small, i.e., 
the /-dependence of K a (l) is similar to one-dimensional case. 
Therefore there is no behavior of spin gap for the total spin 
fluctuation except for very low energy. We note that, for sin- 
gle chain, K a (l) decreases monotonically to K a {l — > 00) — > 1 
and that Ks(l) — 1 for all /. From these /-dependences, it is 
found that a separation of freedoms of charge and spin still 
exists at energy corresponding to formation of the charge gap. 




I 



FIG. 2. The /-dependences of K p (l), K s (l) and K a (l) are 
shown by solid curves, dotted curves and dashed curves for 
g 3 = 0.05 ((1), (4), (7)), psc ((2), (5), (8)) and 0.3 ((3), (6), 
(9)), respectively where parameters are the same as those of 
Fig. 1(a). 
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FIG. 3. The critical values, gsc, are shown as the function 
of £/ep for g = 0, 0.2 and 0.4, where confinement (deconfinc- 
ment) is obtained for §3 > <?3 C (33 < jbc). In the inset, the 
logi/eF-log<73c plot is shown for g = 0.2 (1) and 0.3 (2). 



FIG. 4. The critical values g$ c as the function of 2<?2 — §1 
for t/e F = 0.1 and 0.01 with the fixed 31 = 0.2. The solid 
curves denote the boundaries obtained from Eqs. (2.10)-(2.16) 
while the dotted curves are obtained in terms of the expanded 



In Fig. 3, the £-dependence of gs c , is shown for g\ — 
g2 = g — 0, 0.2 and 0.4 where confinement (deconfinement) 
is obtained for §3 > g^ c (53 < <?3 C ). The £/3 C -dependence of t 
for g = is expressed as 



t/ep ~ ifi exp[-7r/4ff3 



(3.1) 



where ifi ~ 1.2. The intrachain interaction enhances the con- 
fined region. The presence of g leads to a different behavior 
for t with small §3. From the inset which is calculated for 
g = 0.2 (1) and 0.3 (2), it turns out that the (^-dependence 
of t for g 7^ is obtained as 



t/e F ~ K 2 (fcc/Ka) 



1/29 



(3.2) 



where K2 — 0.2, K3 ~ 0.2 and gs c <C 1. The range of g~3 C , in 
which Eq. (3.2) is valid, decreases with decreasing g, e.g. the 
upper bound of gs c is given by 0.15, 0.04 and for g — 0.4, 
0.2 and respectively. Coefficients Ki,Kz and K3 have been 
derived numerically since analytical treatment is very com- 
plicated. We note that the magnitude of §3 C is determined 
mainly by the balance between the charge gap created by the 
umklapp scattering and the energy of interchain hopping as 
is shown later. Coefficients K\ , K2 and K3 have been derived 
numerically since analytical treatment is very complicated. 
We note that the magnitude of g~3 C is determined mainly by 
the balance between the charge gap created by the umklapp 
scattering and the energy of interchain hopping as is shown 
later. As t goes to zero, §3 C reduces to zero and then the 
interchain hopping becomes always relevant in the absence 
of umklapp scattering within the present choice of intrachain 
interaction. 



Now we examine the boundary between confinement and 
deconfinement as the function of g\ and §2 with the fixed 
t. From the calculation with some choices of £/ef = 0.1, we 
find it a good approximation that the boundary with fixed 
(73 depends only on 2g2 — §i- Therefore, g^ c is determined 
essentially as the function of 2§2 — gi, i.e., K p . Such a result 
originates in the fact that G p + t c+(l)-term gives a dominant 
contribution and the effect of g a is neg ligibly small for other 
coupling constants in r.h.s. of Eq. (2.16). In Fig. 4, the quan- 
tity <?3 C as the function of 2<?2 — <?i is shown by the solid curves 
with choices of t/eF = 0.1 and 0.01 where <?i = 0.2. The re- 
gion of §3 > <73 C (§3 < <?3c) corresponds to the confinement 
(deconfinement). The dotted curves is explained in section 
IV. With increasing 2§2 — §1 (i.e., decreasing K p ), the region 
for the confinement is enhanced. 

By use oi response functions for order parameters 



(Eqs. (2.17 )-( 2.19 )), we calculate states at finite temperatures 



where a crossover is shown on the plane of §3 and normalized 
temperature, T/eF, (or energy) in Fig. 5. The dotted curve 
denotes the temperature corresponding to the gap for the to- 
tal charge, A, where A = £Fexp[— l s ] and K p (l g ) = K p /2. 
Around the energy corresponding to the charge gap, the to- 
tal spin shows a behavior similar to one-dimensional case with 
the gapless excitation. The solid curve which is obtained from 
£f exp[— l-t] with t(lt) — 1 denotes the crossover temperature, 
below which the state reveals the property of two-coupled 
chains. Such a temperature becomes lower than the bare in- 
terchain hopping energy, t/e F = OJ^idue to the renormaliza- 
tion by the intrachain interaction£j'E3 The dash-dotted line 
denotes a boundary where the decrease of temperature leads 
to confinement (deconfinement) for §3 > §3 C (§3 < cbc). One 
finds the following four kinds of regions (I) ~ (IV) which are 
separated by these boundaries. The dominant state in region 
(I) is one- dimensional SDW. In region (II), interchain hopping 
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FIG. 5. A crossover is shown on the plane of §3 
and temperature T/eF (or energy) where t/ep = 0.1 and 
gi = §2 — 0.4. The dotted curve and solid curve 
denote A(= epexp[— Z g ]) and epe~ lt respectively, where 
K p (l s ) = Kp/2 and t(l t ) = 1. The dash-dotted line separates 
the region of confinement from that of deconfinement. In the 
inset, the phase diagram of confinement and deconfinement is 
shown on the plane of t/eF and A/ef where the dash-dotted 
(solid) curve denotes the boundary for g — (0.4). 




t/e F 



FIG. 6. The critical value gzc as the function of t/eF 
for g = 0.3. The solid curve corresponds to the calculation in 
Fig. 3 and the dotted curve is obtained by use of the expanded 
K v . The dashed curve is the results for g\ = §2 = <?4 = 0.3. 
In the inset, the phase diagram of confinement and deconfine- 
ment is shown on the plane olt/tF and A/ef for g\ =32 = 0.3 
where the solid (dashed) curve is calculated for 54 = 
(34 = 0.3). 



removes the degeneracy of out-of-phase state and in-phase 
state,LJ but the energy is still high compared with the gaps 
for the transverse fluctuations, which develop just above A 
(dotted curve). In this region, one finds a crossover into the 
out-of-phase SDW state. The present calculation indicates a 
short range correlation for the SCd state in a certain domain 
with small gz{< 0.001) and finite temperatures just above the 
dotted curve. This could retain a traq o that jthe ground state 
for g 3 = is given by the SC stateMHO In region (III), 
the gap of the total charge fluctuation develops. At very low 
temperatures, all the fluctuations become gapful due to rele- 
vant interchain hopping. In this case, the correlation of the 
SCd state as well as the SDW state decays exponentially. We 
note that such a state in the limit of low energy corresponds 
to the "C0S0" phase which rhps been obtained at half- filled 
band by Balents and Fisher.LD In region (IV), the gap of the 
total charge is so large that the interchain hopping becomes 
irrelevant leading to the isolated chainsEj and then absence 
of other gaps. The state in this region hasp* resemblance to 
that of the half-filled one- dimensional chain.tj In the inset, we 
show a phase diagram of confinement and deconfinement on 
the plane of t/eF and A/eF in the limit of absolute zero tem- 
perature, where dash-dotted curve (solid curve) corresponds 
to g = (0.4) in Fig. 3. The region for confinement increases 
by the increase of intrachain interaction, g. The ratio in the 
interval region of 0.01 < t/eF < OA is given by A/t ~ 1.1 
(0.7 < A/t < 0.9) for g = (0.4) £3 These results indicate 
the fact that the deconfinement-confinement transition is de- 
termined essentially by the competition between the charge 
gap and the interchain hopping energy. 



IV. DISCUSSION 



In the present paper, the effect of umklapp scattering on 
two-coupled chains has been examined by calculating the 
boundary for confinement and deconfinement as the func- 
tion of umklapp scattering, interchain hopping and intrachain 
interactions. It has been shown that electrons are confined in 
a single chain when interchain hopping becomes smaller than 
a critical value of the order of the charge gap or the umklapp 
scattering exceeds a threshold. 

We discuss the validity of our calc ulatio n of r enormaliza- 
tion group equations given by Eqs. ( 2.10)-( 2.16). Here we 
calculate these equations by making use of the expansion, 

K^(l) = 1 =F G v (l) H (v = p, a, C and S), where the 

initial conditions are G„(0) — g v /2nvF- Since such a method 
leads to a solution with the SU(2) symmetry, we have done 
following two kinds of evaluations. One of them is shown in 
Fig. 4 by the dotted curve which denotes gs, c as a function 
of intrachain interaction with fixed t/eF = 0.1 and 0.01. A 
good coincidence between solid curve and dotted curve is ob- 
tained for 2t/2 — §1 > 0. The other is shown in Fig. 6 by 
the dotted curve which denotes g 3c as a function of t/tF for 
gi = §2 — g = 0.3. There is a small difference between the 
result without the expansion of K v (solid curve) and that 
with the expansion of K v . Thus, such a calculation may be 
justified within the present choice of parameters. 

Instead of treating two-coupled chains with the split Fermi 
surface, Khveshchenko and Rice applied the bosonization 
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method to two degenerate bands and enlarged the parame- 
ter space of the renormalization group to examine the papi 
hopping for the case of the relevant interchain hoppingJij 
We discuss the effect of the pair-hopping process on the 
present result with the half- fill ed b and . Expressing the 
pair-hopping in terms of Eqs. ( |2.2| )-(2.5), we obtain non- 
linear terms consisting of 6 P +, a +, 0c+ and 9s+, which 
are also found in Eq. (2.9)-|except for a new term given by 
gs cos(\/2^p+) cos(v / 2S (T +)!lj When the interchain hopping is 
larger than the charge gap, the pair-hopping becomes rele- 
vant. In such a region, t(l) and K p (l) of the present paper 
exhibit trajectories similar to those of the reference 15 and 
the differences in magnitudes are small, although the gig-term 
gives rise to a visible enhancement of the spin gap obtained 
by K a (l). It is noted that the interchain hopping is no t ren or- 
malized directly by the gg-term since the r.h.s. of Eq. (2.16) is 
determined only by terms including 6c+ ■ When the interchain 
hopping is smaller than the charge gap, the pair-hopping be- 



comes irrelevant due to confinement, e.g., gs ~ at the energy 
of the charge gap. Therefore it is considered that the effect of 
the pair-hopping is small for the boundary between confine- 
ment and deconfinement. 

Now we examine the effect of the forward scattering with 
the same branch, where the coupling constant is given by 514. 
The Hamiltonian corresponding to g± is expressed as 



Hi 



U, k 1 ,p,cr,i u k 2 ,p,-: 



(4.1) 



which has two kinds effects. One of them appears in K p , K a , 
v p and v a which are written as K p = [{1 — g p /(l + g4,)}/{l + 
5p/(1 + 54)}] 1/2 , K a = [{1-3 CT /(1-S4)}/{1 + 3 CT /(1-34)}] 1/2 , 



M(l + <? 4 ) 2 -ffp] 1/2 and v a =u F [(l-g4) 2 - 
g v /2-KVF- Another is the nonlinear terms given by 



fl 2 ] 1/2 with 



2^2 a 2 J dx {sc+.c- cos (V29 c + - cos V26 C - + gs+,s- cos V29 s + cos V26 S - j 



(4.2) 



where gc+,c- ~ —gs+,s- = <74- Equation (4.2) leads to addi- 
tional terms to renormalization group equations for t(l) and 
Kc(l), which are of the order of o(P) and o(t 4 ) respectively. 
There is also renormalization for Gc+,c- and Gs+.s-- How- 
ever we found that, in the present calculation, the latter effect 
given by Eq. ( |4.2| ) is negligibly small. On the other hand, in 
the former case, there is a noticeable effect of 54, which comes 
from the variation of K p . Note that K p (< 1) increases with 
increasing §4, since the effect of g 4 is equivalent to replacing 
9p by <W(1 +§4) in K p . The increase of K p {< 1) reduces 
the renormalization of umklapp scattering, G p +,c+- Then 
one needs larger §3 to obtain the confinement. An example 
with g = 0.3 is shown in Fig. 6 where <?3 C including the §4- 
term (dashed curve) is compared with that with 34 = (solid 
curve). The quantity gz. c is increased by 34, i.e., the suppres- 
sion of the region of the confinement on the plane of t/eF and 
gs. In the inset, a phase diagram of confinement and decon- 
finement is shown on the plane of i/e F and A/ef for both 
§4 = (solid curve) and §4 7^ (dashed curve) . The fact that 
the critical value of A is reduced by <?4 is understood as fol- 
lows. The magnitude of K p (0) for §4 7^ is larger than that 
for §4 = 0, although the differenc e bet ween these two cases 
is very small as for r.h.s. of Eq. ( |2.10| ). Since larger K p (0) 
leads to slower decrease of K p (l), one obtains smaller gap, A, 
in the presence of $4. Thus, it is found that g4 reduces the 
effects of both intrachain interaction and umklapp scattering. 

Based on the present calculation, we comment on electronic 
states for TMTSF and TMTTF salts, both of which show the 



correlation gap above the SDW state. Such a gap is possi- 
ble for a choice, of T/e F — 10 -2 in Fig. 5 where e F ~ 1Q 3 K 
and T ~ lOif.cJ The fact that the plasma edge perpendicu- 
lar ta chains is present for TMTSF salt (absent for TMTTF 
salt)tl suggests the relevance (irrelevance) of the interchain 
hopping. These behaviors of interchain hopping are found 
when the umklapp scattering satisfies go, ^ g^ c for TMTSF 
salt and g 3 > g 3c for the TMTTF salt. Thus the existence 
of dimerization, which increases §3, is crucial to understand 
property of these salts above the SDW state. 

For the metallic state, the above organic conductors may 
be regarded as a doped system rather half-filling since one- 
particle hopping between chains leads to small deviation from 
commensurate one.cj In this case, the metallic state is ex- 
pected with increasing the doping rate as is shown in a one- 
dimensional caseO Actually a crossover from confinement to 
deconfinement has been obtained in the presence of doping 
even for two-coupled chains.u It will be of interest to study 
such an effect of doping on many coupled chains. 
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APPENDIX A: PHASE REPRESENTATION OF Ti 
AND ORDER PARAMETERS 



By making use of a fc ,p iCTii = \ (-l) l c k . ri . rT . + + c k ri ^.^} (i = 1, 
2) and = L- 1 ' 2 ^^^,^^, Eq. @ is rewrit- 

ten in terms of ^)p >a ,n- The terms for interactions are di- 



vided into two parts, which consist of the scattering be- 
tween electrons with same /i and the scattering between elec- 
trons with opposite fi. Since the former part is treated in 
the way similar to the kinetic energy, we examine the lat- 
ter part which is defined as Hint- By defining Vv.o-.n as 
i>'p,<r,fj,{x) = (1/V27ra) exp (ipkpfj.x + i6p, CT , M ), Hint is written 
as 



Hint = 7 £ / dx 



,f] if] if t I i(an)2-KN a 

91 Vp,a^V-p,^Wp,<r,-fi,V-p,-r,- t t e 



+ 91 i>' v \<y, tl i''-v-<y,v,i''v-o-^'-V,< 



i(ptr)2nN d +in 



+ 91 i! 5i ^- p ,„-^^-^- P ,«T^ ^(JXTW* b +l* 

+ 91 ^'lo^'lp-a-^'v-o -V V'-p^.M e l(P<T)2,rJVd 

+ 01 lb' 1 ^ lb' lb' -i(pcr)27v[N b +nN d ]+in 

+ 01 ^ l/>' - l ( PCT )2 7 r[iV b +p M iV c ]+ l . 

^ i/1 Yp,a,n Y—p,—a^—fx Yp, — a,fj. Y — p : cr, — }i c 

i y2 Yp,a,fj, Y—p,<j 7 —fj. *r — P,cr,fj, Yp,cr,—fj, c 



' i/2 Yp,a,fi Y—p.,—a,—fA Y — p, — cr,/j, YP,cr,—fi c 



T </2 Yp,a,fJ, Y — p,a,n Y — p,a, — fj, Yp,cr,—fj, c 

I i/2 Yp,<7,fi Y — p, — a,fj, Y — p, — a,— fj, Yp,cr,—p, c 

^ i/3 C Yp,a,n Yp, — a,fi Y — p, — <J,[i Y — p,cr,H c 

^ y3 Yp,cr 7 fJ, Yp, — a, — fi Y — p, — (T, — ^ Y — p,a,^i c 



i(afj,)27rN a 



+ no Jp^f^j/I 7 //t // // p ip27r[5 p/i)+ {4Ar a -Ar 7 -Ar 8 }+5p^ _{iV 3 +JV 4 }] 

i/3 c Yp,o-,fi Yp, — a,n Y—p, — o-, — fj, Y—p,o-, — fj, c 

l 



(Al) 



where iV a = [(iVi+7V 2 ) + (iV3+iV4) + (iV 5 +iV 6 ) + (iV7+iV g )]/4, 
JV b = [{{Ni+N 2 ) + (N3+N 4 )}-{(N 5 +N 6 ) + (N 7 + N s )}]/4, 
N c = [{(iV 1 +7V 2 )-(iv3+iv4)}+{(iV5+iV 6 )-(iv 7 +iV 8 )}]/4and 
Na = [{(Ni + N2)-{N 3 + N4)}-{(N5 + N 6 )-(Nt + N 8 )}]/4. 
By defining Nj as the eigen value of Ni, it is found tha t th e 
factor in Eq. (Al) commutes with the Hamiltonian, Eq. (|2.l[ ), 
when N a , Nb, N c and Nd are integers. Such a choice of Hilbert 
space leads to negative sign for 2, 3, 5, 6, 8, 9 and 12-th terms 
in Eq. (Al). By expressing ibp >atlJl (x) in terms of th e phase 
variables, we obtain the non-linear terms in Eq. (2.9). 



Next we examine order parameters. For order param- 
eter, OlsdW|| out , which is expressed as OLSDWy out = 

we evaluate correlation function given by (Olsdwh out ( x ) 
Ol SDW (0)}. By noting that a typical term of this cor- 

I |OUt 

relation function is rewritten as 



(ibi(x)Mx)) (v4(o)^(o)) f 



[fa?{x) fa{x) V?(0) ^(0) e -«-[^+*.]-H- 



= -((^(a!)^(x)) (^(0)^(0))^, 
the correlation function is rewritten as, 

( [- (v4 ^ + ^ Vs) + (v4 + v4 v*)] 

x [- (i>\ fa + fa\ fa) + (V4 fa + fa\ fa)} f 

= {[{faH' 7 + fa2fa) + {fa}fa + ^U'e)] 

x [(faUr+i>Us)+tyU'8+i>u'6)y 



(A2) 



(A3) 



Therefore OlsdWii out in the response function can be ex- 
pressed as 



Olsi 



CT,fX 



(A4) 



which leads to Eq. (2.17) with fa pcr ^(x) = (l/y/2-ira) 
exp (ipkiFiiX + i0 PlCTjM ). The other order parameters are ob- 
tained in a similar way. 



APPENDIX B: DERIVATION OF 
RENORMALIZATION GROUP EQUATIONS 

We evaluate resppssei functions by use of the renormaliza- 
tion grou p method.Ej'E3 By treating the nonlinear terms in 
Eq. ( |2.9[ ) as the perturbation, the response function for 8s+ 
field is calculated up to the third order as 



J 



(T r exp [(i/V2)0 a +{xi,Txj\ exp [-(i/s/2) 9 s+ (x 2 , r 2 )] ) 

: -(K s /2)C/(rf-rF) 



+ 



E l \^ f rf 2 ^ d 2 rl c -(K s /2)U(rf -r*) c -2K„U(r!(-rZ) 

(4n) 2 v 2 J a 2 a 2 

u = p,cr e— ±1 



'I) ^ e «JCs{^('-F-Tf )-l/(rf-ff)-tr(tf-Tf)+t/(r|'-rJ')} _ 



+ (4tt) 2 51 



<: = ±1 ' 



+G 2 + , s _e~ 

' dVf dVf e -(K s /2)C/(rf-rJ) 

a 2 a 2 



-P_ r F ) _ t/ P (T .F_ r P ) _ t/(r |_ r P ) + [/(r P_ r F )} _ ^ 



{ 



(^, s+ e- a ^^-^co6 2»( IS -i*) + ^_, S+ e-'™^- p « P >} 



xe 



-2K s i7(rf 



^ K s{U(rf- r F)-U(rf- r F )- C /( r F- r F) + C /(r-f-rF)} _ A 



{G 2 c+lS _e- 2 ^^'-3 ) CQS 2?o(x3 _ Xa) + cP^e-WKoWtf-rS) j 



xe 



-(2/K s )U{rl-rl) f ^{U (r 



') ^e{C/(rf-rf)-C/(rf-rf)-y(r|'-rf )+C/(rf-rJ)} _ j\ 



_ \^ 4 \^ f d^f d 2 < rf 2 rg c _ ( K R /2)l7(^-r?) c -2i^C/«-r ? ) 

•^-^ (4rr) 3 {;2 ^— ^ / a 2 a 2 a 2 



£ = ±1 ' 



|G„ +J c+G„ + ,s+Gc+,s+e 2KcC/(r 3 r s ) cos 250(2:3 - a? 5 ) 

+G„ + ,c-G, + ,s + G c _, s+ e-( 2 /^)^f-'f)} e - 2 ^^3 F -r) 
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X ^ K s{U(rF -r^-Uir* -r*)-U(rF -r*) + U(r* _ r F )} _ ^ 

-{G„ + ,c + G„ + ,8.G c+ ,s-e- 3KeE ' (r '- r ' ) cos 290(2:3 -z 6 ) 

rt ^ Q -(2/jf c )c/(r!'_ r F)\ 
+Gy+ : c-(j-y+,s-LTC-,s-e 3 5 ' > 



(2/Jf s )[7(rf-rf) 



(•' 



e {C7(rf-rF)-C/(rf-rF)-[/(rf-rF) + C/(r|'-rJ)} 



+ • 



(Bl) 



r 



where U{r" — rj) 



In 



sj (xi - Xj) 2 + vl{n 



for 



i,j = 1,2,3,4, d r v — v v dxdr (y — p, a and F) and 
qo = 2t/v-p. In order to obtain scaling equations of the cou- 
pling constants up to the second order, we need to calculate 



response functions up to the third order for the nonlinear 
terms- Note that these terms do not exist in one-dimensional 
case.Ej By putting rs = r± + r and rs = r^ + r, and expanding 
near r — 0, we obtain the following renormalization in terms 
of effective quantities, 



Kf 



■*sff 2 



-,eff 

r Cp,Sp' 



2 X/ G ^+'S+ K s 



d^_ 

a 



3-2K„-2K s 



dr^ 

a 



p \ 3-2iC„-2/K s 

r 

Q 



dr^ 

Q: 



3-2if c -2X s 



Jo(2q r 



V 2 



dr^ 

Q 



3-2K c -2/K s 



Jo(2q r F ) 



1^2 ^2 



dr^_ 

a 



F s 3-2/K -2K s 



+ 2 G c-,s- 



dr v 
a 



p x 3-2/K -2/K s 



c+G 



dr^ 

Q 



1-2K C 



-2G C - 



.Sp 



G„+,c-GW,sp / — 



Jo(2q r 
FX i-2/jr c 



^ — ' v v I a \ a 



r 



(B2) 



(B3) 
(B4) 



where r" = (x 2 + (v^r) 2 ) 1 ^ 2 , v v = v u /vf and p, p' — ±. The 
second and third terms of r.h.s. in Eqs. ( |B3|) - (|B4|) are ob- 
tained by exponentiating the third order terms-pf Eq. (Bl). 
For the transformation given by a — > a' — ae dl £3 these quan- 
tities are scaled as 

Kf (K',G',q' ,a') =Kf(K,G,q ,a), (B5) 
Gl p y p <{K' ,G' ,q' ,a) = G^,y p i(K, G,q ,a) 

x (a '/a) 7 ""'"'"' , (B6) 

where A", G' and qq denote renormalized quantities. The 
exponent y up y p > is given by y up y p > = 2 - JC* - A'^, . 



By applying this infinitesimal transform to Eqs. (B2), (B5) 
and (|B4|), w e obt ain re norm alization equations given by 
Eqs. (|2.12| ), (|2.14| ) and ( 2.15 ), respectively. In a similar 
way, the renormalization group equation for K V {1) (y — p, 
a and C) is calculated from the response function given by 
( TTe (.i/V2)ev±(.xi,T 1 ) e -(i/V2)0„ ± (x 2 ,T 2 )} and the equations for 

Gv+,c±(l) and G„+,s±(0 [y — p asd a) are calculated from 
the response function for 6„± fieldiij We note that, in case of 
t = 0, thee« equations become equal to the one-dimensional 
equations." 

The renormalization equation for t(l) is derived by evaluat- 
ing the difference of the density between two bands (p = ±), 
which is given by 
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An = 2(/cf+ — fcF_)a + / dxdr ^/cf+ — &f_ / " 



= — 2q a + V2~a J dxdr {d x 9c+{x,t)) 

= -2q a+ 22 —Gv+,c+K c -jr / dxdr (a; sin (V26 c + - 2q x) cos- v / 20„+) 

U — p,<7 

+ / J — Gc+,s P Kcj- / dxdr (x sin (\/20c+ — 2<? cos V^Sp) 
P =± a J 

V — p^G 



Ji(2gor 



> F / F \ 2-2if -2ifP 

dr I r 



Ji(2q r t ) + 



(B7) 



where 4(fcp + — fcp_)/27r = ^ p a /j, tpp^^tpp,^^' The assumption of scaling invariance of Eq. (B7) with respect to infinitesimal 
transformation a' — ae leads to Eq. (2.16). 



The response functions are calculated in terms of the solutions of Eqs. ( p.l0[ )-(2,16). The response function for Oss^ in , 
which is defined by Rss ± ia (xi — x%,Ti — t%) = {T T Oss ± in (^i, r i)Og S± (x%, T2)) , is calculated by writing Rss ± in (x,r) = 

4°L,Jx,t)F SSi . d (x,t) where R^ ± Jx,r) = (a/ ^ + v^) 1/2Kp ' (a/ V* 2 + vW) K ° /2 (a/^/x* + v^f /Kc+Ks)/2 . 
By assuming the scaling relation F B s^(r,a(l),Ku(l),G(l)) = Isa ±M (dl, K v (l), G(l)) ■ F ss± , in (r, a(l + dl), K v (l + dl), G(l + dl)), 
the multiplicative factor Iss± in is written as, 



Iss x>ia = ex P -G a+> c-dl + G a+>s +dl + G C -,s+dl 

+^2 |-G^ + ,c+ M2q a) - G 2 P+ , C - - G 2 P+>S+ - G^+,s-|t/K - A) dl 
+^Kl\Gl + , c +M2q»a) + G 2 +, c _ + G 2 +, s+ + G*+,s- Jtf(r? - r° 2 )dl 
+~{ Yl (-Gl+,c+M2q a) + Gl +tC ./K 2 c + G 2 V+S+ K\ - G 2 V+ , S - J 
+G 2 c+! s+ (-1 + K 2 S ) M2q a) + G 2 C+ ^ (-2) J (2q a) 
+G 2 c_, s+ (l/Kl + Kl) + G^_, s _ (l/K 2 c ~ 1) j U(rf - r 2 F ) dl 



which leads to -Fss^ in expressed as 

Fss ±in (r,K, G) = exp 



'ln(r / ol) 

£ hx[l ss±M (dl,K(l),G(l))] 



(B8) 



(B9) 



We note that terms including the second order of the coupling constants are rewritten in a simple form. For example, one 
obtains 

\J dl {j2 (-Gl + , c+ J (2q a) + Gl +tC _/K 2 c + G 2 V+ , S+ K 2 S - G'+.s^ + G 2 C+ , S+ (-1 + ^ s 2 ) J (2q o a) 



^2 . r /2\ , ^2 





r 




[a(l)\ 
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dl 



dK c {l) dK s (l) 



K%(1) dl 



dl 



In 



a(l) 



K c (0) 



+ ifs(0)}la(l) 



ln(r/a) . f 



+ K s (l) 



(BIO) 



where = ae l . Thus, the normalized response function, iiss^ in ( a; ) r ) (= -Rss^ in ( a; >' r ) ' 2(7r«) 2 ), is expressed as, 



^ssj. in {x,r) = exp 



In yj x 2 + (v p t) 2 /a 



dl 



2K p (l) 



x exp 



x exp 



fin \J x 2 + (v (t t) 2 I ol 



/ dl \ --K„{1) - — G« r+ ,a-(/) + — G„+, s +(0 f 

J I 2 u CT w CT J 

In -y/ ic 2 + (^f t ) 2 / a 



dli-- 



K c (l) 



+ K s (l) + G C -,s+(0 



(Bll) 



which leads to Eq. (2.22). Other response functions are obtained in a similar way. In deriving Eqs. ( 2.20 )-( 2.2S ) , we replaced 
v p and Va by vf, which may cause a slight deviation of the numerical factor. 
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